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Abstract 



We show that the logarithmically rising static potential between opposite-charged 
sources in two dimensions is screened by dynamical fields even if the probe charges are 
fractional, in units of the charge of the dynamical fields. The effect is due to quantum 
mechanics: the wave functions of the screening charges are superpositions of two bumps 
localized both near the opposite- and the same-charge sources, so that each of them 
gets exactly screened. 



1 Introduction 



The static potential between trial external charges, or the Wilson loop expectation value, 
carries important information about infrared behavior of gauge theories. An infinite growth 
of the potential provides the simplest criterium for confinement. However, in certain theories 
the rising potential can be screened at large distances by dynamical fields. The screening 
is inevitable if dynamical charges can form neutral bound states with external sources, for 
example, when the external and the dynamical charges belong to the same representation of 
the gauge group (have the same magnitude in the Abelian case). Since the potential between 
bare charges can be arbitrary large, at some point the creation of a pair from the vacuum 
becomes energetically favorable. Each of the created charges couples to the static charge 
of the opposite sign. The interaction energy between resulting bound states no more grows 
with the separation. 

The question of whether fractional charges can be screened or not is more involved. This 
problem has been studied in two dimensions, both in Abelian [p], [2|, [3| and in non- Abelian 
[0. |J models. It appears that massless matter fields can screen any Abelian fractional charge 
[p], §. In non- Abelian case, massless fields in any representation of the gauge group screen 
sources in the fundamental representation 0, §], as follows from comparison of the models 
with massless adjoint matter and with multiple flavors of fundamental matter ||. 

In this paper we consider the problem of charge screening in three-dimensional scalar 
QED. This theory is confining when matter decouples, as the Coulomb potential in two di- 
mensions grows logarithmically with distance. We consider bosonic theory to purify the dis- 
cussion, because fermions, at least massive, screen any charge in 2+1 dimensions ||. Three- 
dimensional fermions induce the topological mass for a photon at one loop thus changing 
the logarithmically rising Coulomb potential to the exponentially decreasing Yukawa one. 
This phenomenon does not take place in scalar QED, in which the photon remains massless. 

We argue that, nevertheless, any fractional static charge in 3D scalar QED is screened, 
at least in the weak coupling regime. The coupling constant e 2 in three dimensions has 
the dimension of mass, so the weak coupling means that the ratio e 2 /m is small. We do 
not expect any abrupt changes to happen as this parameter is increased. Therefore, the 
screening, most probably, persists in the strongly coupled theory as well. 

We consider the static potential between external charges of the magnitude q separated 
by a distance L (we imply that < q < 1 for simplicity - the screening of the integral 
part of the charge is obvious). The screening mechanism is very simple and is based on 
the consideration of a two-particle state. If g 2 e 2 ln(L/i) > 2m, the interaction energy is 
sufficiently large to create a pair of dynamical charges from the vacuum. The parameter [i 
is determined by the typical scale of the screening charge distribution, actually fi ~ \J me 2 . 
At first sight, such rearrangement of the vacuum can only change q to q — 1, but cannot 
stop the logarithmic growth of the potential with L, since the dynamical and the external 
charges form a bound state which is also charged. However, the wave functions of dynamical 
particles need not be localized near the opposite-charged sources only. Suppose the wave 
function has the form of a superposition of two states well localized near each of the static 
sources. Let these bumps be normalized to the probabilities p and 1 — p, respectively. 

This situation is schematically illustrated in fig. |l|. It is clear that if p — (1 — p) = q, that 
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Figure 1: The delta-peaked external charges (dashed lines) together with the distribution of 
screening charges (solid lines). 



is, if p = (1 + q)/2, the net charge localized near each of the external sources sums up to 
zero. At the same time the total probabilities to find the positive- and the negative-charged 
particles are both unities. Only short-range interactions are present in such configuration of 
charges; its energy does not grow with L, and, thus, it becomes energetically favorable at 
very large separations. We will argue that leaking of the charge to the region where it is 
classically repulsed actually takes place for the Klein-Gordon particle in the electric field of 
well-separated static charges. After that we show that the configuration described above has 
smaller energy than bare sources for sufficiently large L. The stability of this configuration 
can be heuristically explained as follows. The logarithmic Coulomb potential of the point 
charge is singular at short distances, so it is quite natural that the negatively charged particle 
form a bound state with the positive external source. This bound state is described by a 
larger part of the double-bump wave function of type plotted in fig. [I]. The charge of this 
bound state is —(1 — q) /2, SO , cLS £L whole, it attracts the positively charged particle, which 
explains the stability of the smaller part of the wave function. 

The paper is organized as follows. In Sec. 2 we diagonalize the Hamiltonian of the scalar 
QED in the presence of external charged sources in the weak coupling approximation. In 
Sec. 3 we discuss the Klein-Gordon equation for a charged particle in the electric field of a 
dipole. In Sec. 4 the two-particle state becoming energetically favorable at large separation 
is considered in more detail. In Sec. 5 and in Appendix we comment on the path integral 
for scalar QED in the presence of external charges. 
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2 External charges in scalar QED 

We consider scalar QED in three dimensions. The Lagrangian density of this theory is 

C = -- F^F^ + {D^D^ - m 2 &<5>, (2.1) 
where $ is a complex scalar field and 

D/x — dp + ieA„. (2.2) 

For our purposes the canonical formalism is more appropriate. 

In the Schrodinger representation, A = and E { = F 0i = A u A i} II = $, W = $, 
are canonical variables: 

[A(x), J E? j (s/)] = iM(^-2/)> ( 2 - 3 ) 

[<£>(*), n(y)] = iS(x -y) = [&(x), nt(y)]. (2.4) 

The Hamiltonian is 



H = J d2x l E i + \ F i + 11+11 + (A^) f A^ + ^ 2 ^ f $ 



(2.5) 



The physical states in the presence of external charges are subject to the Gauss' law con- 
straint: 

<9^|^phys> = (J + P)|*phys), (2.6) 

where Jo is the charge density operator, 

J Q (x) = ie(&(x)ltf(x) - $(x)Il(xj) (2.7) 

and p is the density of external sources. In our case of two well-separated point-like charges 

p(x) = qe{5 (x + L/2) - 5 (x - L/2)) . (2.8) 

The potential of interaction between charges is equal to the difference of the ground state 
energies in the sectors of the Hilbert space defined by the Gauss' law with and without 
external sources. 

The qualitative picture of charge screening is based on purely classical notion of electric 
field which becomes strong enough to create pairs. It is difficult to visualize this picture 
in the Hamiltonian formalism, where electric fields are operators in the Hilbert space and 
logarithmically rising electrostatic potentials are somehow encoded in the dependence of the 
wave functional on A t . However, it is possible to introduce classical, c- number fields which 
play the role of electric potentials in the Hamiltonian formalism, despite A = by definition. 

Consider the following Hamiltonian: 

H(ip) = H + j d 2 x if (diEi -J -p), (2.9) 
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which now acts in the unconstraint Hilbert space. The eigenfunctions and the eigenvalues 
E of this Hamiltonian are the functionals of ip: 

H((p)\V) = E\V). (2.10) 

We want to show that, if (p is determined by stationarity condition 

8E 

the state |\&) satisfies the Gauss' law (|2.6|) and is an eigenstate of the Hamiltonian Q2.5Q with 
the eigenvalue E. 

Both of the Hamiltonians ( |2.5|) and ( |2.9|) commute with the Gauss' law and with one 
another, so they can be simultaneously diagonalized. Therefore, it is sufficient to show 
that the Gauss' law is satisfied in average. This immediately follows from the stationarity 
condition (|2.11|) and the Schrodinger equation fl2.10|) : 

6p (#|#) ($|$) ($|$) 

The interaction of the scalar fields with photons entering the Hamiltonian through the 
covariant derivative squared term can be disregarded in the weak coupling limit since it is 
of order e 2 /m and is not enhanced by a ln(L/x) factor. The Hamiltonian H(ip) is quadratic 
in this approximation: 

1 „o _ „ 1 



H{tp) = J d 2 x 

+n f n + iep($Il - $ f n f ) + d^d^ + m 2 $ f $ - pp] , (2.12) 



- Ef - Eidip + - El 



and can be explicitly diagonalized. 

Solutions of the Schrodinger equation for the Hamiltonian Q2.12 ) have factorized form 

matter 

[$,$ t ]. (2.13) 

We first consider the gauge-field part of the wave function. The ground state is described 
by a Gaussian wave functional: 

^gau ge [^] = exp (-- J d 2 xd 2 y A i (x)K ij (x,y)A j (y) + i J d 2 x £ l {x)A i (x)^j . (2.14) 

Substituting this expression in the Schrodinger equation (electric fields act on the wave 
functional as variational derivatives: E± = —i5/5Ai), we obtain for and Sf. 

djdj 

-d 2 



K iS = [6ij + ] , (2.15) 



Si = ditp. (2.16) 
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The energy of this state is 



(2.17) 



where E is the divergent zero-point energy 



(i 2 p 



•P, 



which does not depend on if and is omitted below. 

The Hamiltonian for matter fields can be diagonalized introducing creation and annihi- 
lation operators: 

[#(¥>), a + ] = £a f . (2.18) 
Since Hamiltonian commutes with electric charge 

Q = J d 2 xJ , (2.19) 

operators satisfying eq. ( 2.18Q are linear combinations of II and $t ( or Q f ]j\ anc [ <j>) : 

d 2 x (llif) + $ f ^) . (2.20) 



Substituting this operator in eq. ( 2.18D we find that if) = i(e + eip)if>, where if) and e are 
determined by the equation 



(e ±n + e^f + d 2 - m 2 \ if) ±n = 0. 
Here the subscripts ± mark positive- and negative-energy states: 

e +n > 0, £_„ < 0. 



(2.21) 



(2.22) 



The equality ( |2.21p is nothing but the Klein-Gordon equation for eigenmodes in the time- 
independent external field = S^Qtp. Its solutions form two complete sets of functions 
normalized by l| 



J d 2 xif)± m (e± m + e± n + 2&p) ip± n = ±5 mn , 



d 2 x if)± m (e ±m + e Tn + 2eip) ip 



-f n 



0. 



(2.23) 
(2.24) 



a+, 
at, 



These eigenfunctions determine two sets of operators 

- 1 d 2 x U(x)ip +n (x)+i^(x)(e +n + e(p(x))ip +n (x) 
d 2 x if) +n {x)U\x) - iip +n (x) (e +n + eip(x)) $(x) 
ip- n (x)U ] (x) - iif)-„(x) (e^ n + e(p(x)) $(x) 
U(x)if)- n (x) + i&(x) {e- n + eip(x))if)- n (x) 



a_ 



d 2 x 
d 2 x 



(2.25) 
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which create and annihilate particles of charge ±e and energy |e± n |: 

[<2,(4„] = ±ea^ n , [Q,a ±n ] = Tea± n (2.26) 
[H(<p),al. n ] = ±e± n a ] ± n = \e± n \a}± n , [H((p),a± n ] = =R:± n a ±n = -\e± n \a± n , (2.27) 

[a± n , a±m] = finm, [d± n, d± m ] = = [a± n , % m] • (2.28) 

The field variables are expressed in terms of creation and annihilation operators as 
$(x) = (i) +n {x)a+n - ip- n (x)aL n ) , 

n 

n(^) = XI [ a +n^+n(^) ( £ +n + e<p(x)) - a_ n ip_ n (x) (e_ n + e<p(x)) . (2.29) 

n 

The total energy is comprised of -E gail ge given by eq. flUTD , the energy of the matter 
fields, ^matter, and the source term: 



E = - d 2 x 



- (^) 2 + 



+ -^matter- (2.30) 



The energy of the state satisfying the Gauss' law corresponds to the extremum of this 
functional, which is determined by the following equation: 

-cP<p + p+ (J ) =0, (2.31) 

where we used the fact that SE^^/Sip = (5H matter ((p)/5(p) = — (Jo). Note that this 
extremum is a maximum of E(ip). 

When the separation of external charges is not very large, an empty state, 

a±„|0) = 0, (2.32) 

has the lowest energy. The vacuum charge density can be expanded in powers of 1/m: 

e 2 

(OlJolO) = const — d 2 (p + ..., (2.33) 
m 

and is small compared to the first term in eq. ( |2.31|) . Therefore, we can safely omit vacuum 
contributions to the energy and to the charge density. 
The equation ( |2.31| ) is then solved by 

eq \x + L/2\ 

^T, hi \^m (2 ' 34) 

and Q2.30D gives the Coulomb law: 

n 2 p 2 T 

Here ( is an UV cutoff necessary to regularize an infinite Coulomb self-energy of the static 
charges (see sec. 4 for the precise definition). However, the logarithmic raise of the interaction 
energy cannot last infinitely without a substantial rearrangement of the vacuum. To study 
it, we consider next the spectrum of the Klein-Gordon equation for the potential fl2.34| ) of 
the electric dipole. 
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Figure 2: The potential ( ggg ). 
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Figure 3: Spectrum of the Klein-Gordon equation for the potential (|2.34 ). 



3 Klein-Gordon equation in the dipole field 



First, when the distance L is small, the Klein-Gordon equation ( |2.21 ) has no normalizable 



solutions and its spectrum consists of two continua with e + > m and with < —m. 
Near the boundaries of the spectrum, for e± = ±(m + A), A <C m, the non-relativistic 
approximation can be used: the first term in eq. ( |2.21|) can be expanded in ey?, and the 
eigenfunctions ip± satisfy ordinary Schrodinger equation with the potential ^feip: 

^-9 2 T e^-A ± )^ ± =0. (3.1) 

The potential energy for positively charged (coming from the upper continuum) particles, 
shown in fig. ^|, has a form of the separated peak and the well. For sufficiently large L the 
attraction by the well becomes strong enough for a discrete level to appear. Since a positive 
charge is attracted to the source at x = L/2 and is repulsed from the one at x = —L/2, its 
wave function is localized near x = L/2. The wave function of the negative charge is 
localized near x = —L/2 and its energy is e_ = — £+o by symmetry. The lowest positive- 
and negative-energy levels converge with the increase of L and collide at zero for some critical 
value of L = L (fig. |3|). After that they do not disappear, but rather go off to the complex 
plane §|. This behavior of the eigenvalues is generic for Klein-Gordon equation in strong 
electric fields [§, |9|, [UJ. For L > L vacuum polarization can no longer be neglected and 
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the external electric field creates a pair of charged particles in the vacuum. This effect is 
analogous to a pair creation in the field of a heavy ion with the nuclear charge Z > 137 0. 

For large L the ground state wave function is no longer localized near the well of the 
potential. Rather, the wave function has the shape with two bumps, like the one used in the 
qualitative arguments in the introduction. As L grows, the charge leaks from the well to the 
region where the potential is peaked. This, at first sight, anti- intuitive behavior reflects a 
generic property of a charged Klein-Gordon particle to form a bound state in a sufficiently 
strong repulsive electric potential [|, ||, [Hj . The fact that the charge is redistributed between 
the well and the peak of the potential can be proved by the following arguments. For L = L 
both eigenvalues e_ and e +0 are equal to zero and the Klein-Gordon equation formally has 
the form of the Schrodinger one: 



-d 2 




x + L/2\Y 2 
1 + m 



x-L/2\J 



±0 



(L 



(3.2) 



The potential here has the shape of a symmetric double well. The ground state wave function 
is symmetrically distributed between the two wells. By continuity reasons, as L is decreased, 
the charge begin to leak from the region near the repulsive source to the attractive one, and 
eventually all the positive charge is concentrated near +L/2 and the negative one near —L/2. 



4 Screening of the logarithmic potential 

So far, the vacuum sector was considered. However, when e 2 ln(L/z) ~ m, the two-particle 
state 

|2> = 4 aL |0) (4.1) 

can become energetically more favorable, if the dynamical charges screen the sources and 
reduce the Coulomb energy by an amount sufficient to create a pair. The energy of the state 

-^matter = £+0 — £-0 (4-2) 

and the induced charge density is 

(2| J \2) = 2ei> +0 (e +0 + e(p)ijj + o + 2e^_ (£-o + e(p)ip- - (4.3) 

Vacuum contributions are neglected here. 

The charged particles which cause the screening of the static charges are non-relativistic, 
since the screened electric fields are small everywhere, unlike the unscreened ones. Therefore, 
it is possible to use the non-relativistic approximation (|3.1|) to the Klein-Gordon equation. 
In this approximation, the wave functions tp± - we omit the subscript for brevity - are 
normalized to l/(2m), as follows from equation (|2.23|) . For the sake of clarity it is, however, 
convenient to introduce new wave functions, ip' ± = V 2mip± , normalized to unity, as it is 
custom in the non-relativistic limit: 

/dW±| a = l, (4.4) 
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whereas the induced charged density is 

(Jo) = e(K| 2 -K-i 2 ) 
The total energy of the two-particle state is 



(4.5) 



E = 2m + — fd 2 x(\d<4)' + \ 2 + \d<4)'_ 



1 

47T 



d 2 xd 2 y \p(x) + (J (x))] In \x-y\ \p{y) + (J (y))} ■ 

(4.6) 

in eq. (|2.30| ), and 



This expression is obtained after 2m + A + + A_ is substituted for E m[LtteT 
the solution of the Poisson equation fl2.31| ) is substituted for <p. 

The energy (|4.6|) can be regarded as a functional of ip' ± . It is straightforward to check 
that the minimum of this functional is determined exactly by the Schrodinger equation ( |3.1| ). 
The coupled set of equations (|2.31|) , (|3.1|), (|4.4j) and ( |4.5|) constitute a rather complicated 
eigenvalue problem. The ground state corresponds to the global minimum, and can be found 
numerically. Instead of solving these equations directly we will suggest simple variational 
wave functions ip'± corresponding to a two-particle state whose energy does not grow with 
the separation of the external charges L. 

We take the variational wave functions in a form of a superposition of states with charges 
±e localized in the vicinity of the sources with charges ±eg, q < 1, see fig. |I| For symmetry 
reasons we take ip'_(x) = ip' + (—x): 



ip' + {x) = 0i (x - 

ipU x ) = 02 (x - ^) + 01 ( x + f ) ■ 



The functions 0i,2 are supposed to be well-localized and normalized to 



d 2 x 4> 2 



1=F9 



(4.7) 



For large L the overlap between 0i and 02 is exponentially small and can be neglected, 
therefore the wave functions ( f4.7|) are normalized to unity. This ansatz corresponds to the 
distribution of charges described in the introduction. The (l + q)/2 portion of the dynamical 
charge is localized near the source of the opposite sign and the (1 — q)/2 portion is localized 
near the one with the same sign. 

Neglecting the exponentially small overlap of 0i 2 we get for the total charge density: 



p(x) + (J {x)) 



(i ,) [x + - 



+ 



L 

X+ 2 



L 

X+ 2 



q8 



x 



L 



x — 



L 



x — 



(4.9) 



2 J ' TL V 2 J rz V~ 2, 
We see that the total charge density noticeably differs from zero only in the vicinity of the 
points x = —L/2 and x = +L/2. The screening of the delta-peaked external sources is 
achieved if the integral of the total charge density over the region much smaller than L is 



zero. This is guaranteed by our choice of the normalization condition ( |4.8|) . 



10 



To make an estimate of the minimal energy of the two-particle screening state and, hence, 
of the critical separation between the external charges where the rising potential breaks up, 
we take a simple Gaussian ansatz for the wave functions </>i 2' 



<t>iA x ) 



1=F9 



exp 



x 



4a?, 



(4.10) 



where the widths of the wave functions a 12 are the variational parameters. To make all 
integrals finite we shall temporarily introduce a Gaussian smearing of the external sources 
replacing 



qS(x) 



q 



exp 



x 



C^o. 



(4.11) 



2ttC 2 " V 2CV 

Now one has to substitute the trial wave functions ( [4.10| ) into the energy functional Q4.6| ) 
and to find the best widths a^ from its minimum. The integrals are readily performed by 
using the Fourier transforms. Recalling that the Fourier transform of ln(x 2 )/47r is —1/k 2 
we get for the interaction or the potential energy term in the total energy (the last term in 
eq.flO)): 



Epot(L) 



d 2 k e ik ' L 



(27T) 



k 2 



qe 



2 fc 2 l-o 4 k2 
— _| —L e — 



1 + q 4 k2 
— - — e 2 



(4.12) 



The term proportional to exp(ik ■ L) accounts for the interaction between the regions near 
x + L/2 and near x — L/2, while the subtracted term proportional to unity takes into account 
the self-interaction of charge distributions inside these regions. 

Eq. (|4.12 ) should be compared to the interaction of two bare external charges: 



Ebare(L) = — e z 



d 2 k e ik ' L - 1 
(2?r) 2 k 2 



qe 



e 2 q 2 

47T 



' L 2 s 



(4.13) 



Terms exponentially small in L/C, have been neglected here. The integral (^4.12|) is immedi- 
ately calculated using ( |4.13|) , yielding 



E. 



pot 



e 

Att 



L 2 



l-g^ 2 L 2 

ln 4o| + 



l + q\ 2 , L 2 

4a| 



+2q- 



q 



In 



2q 



l + q 



In 



L 2 



2^i±^ln. 



(4.14) 



2 _ "2(C 2 + af) * 2 2(C 2 + a|) 2 2 2(a 2 + a|)_ 

The coefficient in front of In L 2 is zero, so that the energy is now independent of L, up to 
exponentially small corrections which are neglected. It is exactly the screening effect we 
are after, and it is due to the choice of the normalization of charge distributions, eq.( |4.8| ). 
Neglecting also the spread of the external charges ( 2 as compared to a\ 2 we get 



E, 



pot 



f_l r 
47r4 



2(1 - q)(l + q) ln(a 2 + a 2 ) - (1 - q){\ + 3q) lna 2 - (1 + q)(l - 3q) In a 2 
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-2(1 + 3q 2 ) In 2 - 4g 2 In ( 2 . (4.15) 
The kinetic energy term (the second term in eq.( |4.6|) ) is 



1 ( 1-q l+q 

The sum, E kin + E pot , has a minimum at 



^=^Ur + V) (416) 



2 _ 4tt l±g + 4g 2 + (l±g) v ^+87 
ai ' 2 "me 2 4g 2 (l T g) ' 10 

which should be substituted into ( |4. and ( |4.15| ) to get a variational estimate of the energy 
of the two-particle state screening the external charges. Naturally, the energy- at-rest, 2m, 
should be added, too. 

It follows from ( |4.17| ) that the distribution of the dynamical charge having the same sign 
as the external charge (ai) is broader than that having the opposite charge (02). For example, 
if the external charge is one half of the dynamical charge (q = 1/2) the same-charge cloud is 
about 2.5 times broader than the opposite-charge cloud. See the table, where examples for 
other values of q are given. 

Finally, the energy of the two-particle ground state can be written as 

e 2„2 / 4tyC 2 \ 

B = 2m+ iH ta ;^H' (418) 

where C q is a number of the order of unity coming from substituting the best values of 
ai >2 given by eq. ([4. 17|) into E kin + E pot . Notice that the dependence on the spread of the 



delta-peaked external charges, (, is the same as in the case of the bare charges, eq. (|4.13|) . 
This is because the extended dynamical charge distribution cannot screen the logarithmic 
potential at small separations. 



Table 



Q 


_ /me 2 

a iy — 




c q 


0.1 


7.96 


6.53 


0.536 


1 


3.49 


1.85 


0.589 


-f- 


3.19 


1.26 


0.649 




3.44 


0.976 


0.717 


0.9 


5.65 


0.776 


0.824 



The logarithmically rising potential between external charges at large separations breaks 
up when the bare energy ( 4.13 ) exceeds the energy of the screening state, eq. flngp . It 
happens at the critical separation between the external charges 



Cn 



I 4tt 
me 2 q 2 



exp 



/ 47rm 

V e V , 



(4.19) 
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Since we assume the non-relativistic limit, m 3> e 2 /4-7r, this distance is exponentially 
large. Notice that the widths of the screening distributions ai 2 as given by eq. (|4.17|) are 
much less than the critical distance L c , which justifies neglecting of the overlaps between the 
screening clouds belonging to the two centers. 

The numerical values of the coefficient C q are given in the table for certain values of q, 

together with the values of the widths a\ t 2 measured in natural units of ^me 2 /An. Since we 
have used a variational estimate for the ground state energy, the true minimum can be only 
lower, that is to say that the numerical coefficient C q of the order of unity in eqs.( |4.18| ),( [4.19| ) 
can be somewhat smaller than given in the table. However, the dependence on the algebraic 
parameters in these equations follow from the dimension analysis, and is of a general nature. 

At L Pd L c the logarithmic growth of the potential stops; more precisely it slowly grows 
approaching its asymptotic value at infinity ( 4.18Q , the deviation corresponding to the resid- 
ual forces between neutral charge clouds at x = ±L. 



5 Path-integral approach 

The arguments we used above are in essence the variational ones. For this reason, we 
preferred to use the Hamiltonian formalism. Although the discussion of the screening from 
the path-integral point of view is beyond the scope of the present paper, we would like to 
outline how the main ingredients of our analysis can be derived from the path integral. Here 
we consider only the vacuum sector. 

The vacuum average of the Wilson loop infinitely stretched in the time direction, which 
determines the energy of two static charges, is given by the path integral 



(5.1) 



Z = J DAD&D& exp ji J d 4 x -- F^F^ + (D^D^<t> - m 2 &<5> 

+iqe J dx° (A (x°, -L/2) - A (x°, +L/2)) } . 

Integration over the scalar fields induces an effective action for the gauge potentials: 

AS = iTr \n(-D 2 - m 2 +i0). (5.2) 

The next step, justified by the smallness of the coupling, is to calculate the remaining integral 
over An in the saddle-point approximation. This amounts to solving classical equations of 
motion taking into account the induced action ( |5.2j ) and the source term in (571) . We are 
going to show that these saddle-point equations are nothing but the ones derived in Sees. |2], |3| 
in the Hamiltonian formalism. 

Owing to the symmetries of the problem the classical fields are time-independent, and 
Af = 0. For Aq we get the Poisson equation, 



where p is the same as in 
5 



-didiJ$ + p+(J Q ) = Q, 
and the induced charge density is 



(5.3) 
(5.4) 



x=y 
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where 



G(x,y) 



x 



—D 2 -m 2 + i0 



y 



(5.5) 



We leave the calculation of the induced charge density to Appendix, where we show that 

ein- Gordon equation ( 
saddle-point equation (|5.3|) coincides with (|2.31|) . 



it reduces to solving the Klein-Gordon equation ([2.2 1| ) with (p replaced by v4q and that the 



6 Conclusions 

To summarize, the infinitely rising potential between fractionally charged external (probe) 
sources is screened in (2+l)-dimensional scalar QED. Of course, if the mass of the dynamical 
fields is large, the rising potential persists at intermediate scales. The critical distance is 
exponentially large in m/e 2 , in contrast to 3D spinor QED where the screening length is of 
the order of 1 /e 2 || . 

The screening is a typical quantum- mechanical effect: the wave functions of the screening 
particles are superpositions of two distinctive bumps localized near external sources of both 
signs and carrying fractional charge. In case of a half- integer charge of the probe (q = \) the 
bumps carry charges | and ~, so that the total probability is unity but the external sources 
are completely screened, see fig. [IJ. 

It is interesting that the screening effect would be probably not easy to observe from 
Euclidean lattice simulations of the theory (given by the partition function ( |5.1| )). Indeed, 
the essence of the mechanism is a formation of two bumps in the screening wave functions, 
which is a kind of tunneling effect. The larger the separation between sources L, the longer 
computer time one would need for this effect to come into action, with the time growing 
exponentially with L. 

Nevertheless, it would be very useful to check the screening of fractional charges by lattice 
simulations, in view of apparent analogies with a more difficult case of non-Abelian gauge 
theories. 



Acknowledgments 

D.D. would like to thank Victor Petrov for many discussions in the past that stimulated this 
investigation. K.Z. is grateful to NORDITA for hospitality while this work was in progress. 
The work of K.Z. was supported in part by NATO Science Fellowship, CRDF grant 96-RP1- 
253, INTAS grant 96-0524, RFFI grant 97-02-17927 and grant 96-15-96455 for the promotion 
of scientific schools. 



14 



Appendix A Eigenfunction expansion of induced 

charge density 

Since the background field is time-independent, the eigenmodes of the operator — D 2 — m 2 , 

(-D 2 - m 2 )u = ku, (A.l) 

have the form 

u=e i£X °u, (A.2) 
where u depends only on spatial coordinates and satisfies the equation 

{e + eA$) 2 + d l d l - m 2 ] u n = n n u n . (A.3) 

The eigenvalues k„ are the functions of e. The eigenfunctions u n are supposed to be nor- 
malized to unity for any given e. 

The Green function ( |5.5|) can be expanded in the eigenfunctions as 



G{ ) = [ + °° ^ e i^-y°) V ^KM, (A . 4) 
v ,yj J-oo 2vr V K n (e)+iO v ; 

The integral over e can be calculated closing the contour of integration in the upper or in 
the lower complex half-plane, depending on the sign of x° — y°. The charge density is defined 
in the limit x° — > y° and, in principle, depends on how this limit is approached. The correct 
prescription is to take the limit symmetrically in x° — y°. 

The integrand in ( A.4 ) has a pole if K n (e) = 0. The values of e which satisfy this condition 
are determined exactly by equation ( |2.21| ). The only difference is in the normalization of the 
functions u n and ^ n 0: 

V>n = (A.5) 

\NJ 



N n = 2j d 3 x u n {e n + eA c J)u n . (A.6) 
For sufficiently weak fields all poles lie on the real axis. Near a given pole one has 

K n {e) = N n {e -e n ) + ..., (A. 7) 

since, as it follows from eq. ( |A.3| ), 



de 

For weak fields N n is positive for positive e and negative for negative ones. Therefore, the 
contour of integration in eq. (|A.4J) passes the poles on the positive semiaxis from above and 
the ones on the negative semiaxis from below (fig. |]). By continuity, this rule remains valid 
for strong fields until the positive and the negative poles collide. It can be shown that exactly 



*We use here slightly different notations than in the main text where we discriminated between the 
negative- and the positive-frequency modes. 
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Figure 4: The contour of integration in the complex e plane. 



when two poles collide (fig. |3|) and go off to the complex plane the integral iV n turns to zero 

!■ 

As a consequence of eqs. ( |A.7p , (|A.5|) , the functions u n are replaced by ij) n in the residues 
of the integral ( |A.4|) . Taking into account that the covariant derivative Dq acts on the 
eigenf unctions (|A.2|) of the Klein-Gordon operator as i(e + eA^), we find for the induced 
charge density: 

{Jo) = e^1) n (e n + eA$)il> n . (A.8) 

n 

It can be checked that this expression coincides with the vacuum average of the charge 
density operator defined in sec. ||. 

The above consideration implies that all eigenvalues of the Klein-Gordon equation lie 
on the real axis, which is true for sufficiently weak fields, that is, for sufficiently small 
separation between external charges L. For very strong fields (at L > Lq) some poles move 
to the complex plane, as discussed in Sec. 0, which invalidates the analysis above. This 
signals that for large L the ground state is rearranged and the true vacuum is different from 
an empty state with no dynamical charges. 
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